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Abstract 



New features are described for models with multi-particle area-dependent 
potentials, in any number of dimensions. The corresponding many-body 
field theories are investigated for classical configurations. Some explicit 
^ , solutions are given, and some conjectures are made about chaos in such field 

theories. 

in 
o 
t> 
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. Area-dependent potentials, V (A = r\ A r2), or their mathematical equivalents, 

appear in several physical problems of contemporary interest. Notable among 
these problems are Yang-Mills theory (especially for spatially homogeneous config- 
urations, with only time dependence) ||, extended supersymmetric field theories 
with (pseudo) scalar self-interactions [|18|, [13], p20| , and more recently, membrane 
models JTl], [llj. In the last century, Feynman |14[ even assigned to students an 



exercise involving such potentials: Show the energy spectrum for a quantized 
areal potential model is discrete and quantum particles cannot escape from such 
a potential, even though the spectrum for the classical model is continuous and 
classical particles can escape along special trajectories for which A = 0. Many of 
us pondered this problem over the intervening years, especially in the field theory 
context Q . Meanwhile, Barry Simon worked out five or six solutions to Feynman's 
exercise and published them [^5] . 



Classically, such potential models are interesting insofar as they may provide 



simple examples of chaotic systems [f24| , |2j. While there are special trajectories for 
which the motion is quite regular, including those for which particles can escape 
(such as straight line, free particle motion), for most trajectories this is not the 
case. The current consensus is such models are not integrable in the Liouville 
sense, and do not admit the construction of a Lax pair. Is this really so? 

*Bascd on a talk given at the "Coral Gables" conference, Fort Lauderdale, Florida, 12 Decem- 



ber 2001. © American Institute of Physics, http://proceedings.aip.Org/proceedings/#624: Coral 



Gables Conference on Cosmology and Elementary Particle Physics; B. N. Kursunoglu, Editor. 
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Quantum mechanically, there are obviously more interesting questions to ask 
about such potential models than simply whether the energy spectrum is discrete. 
In particular, just what is that discrete spectrum? Is the quantized model, or 
some simple variant of it, completely integrable [O, even if the classical is not? 



Do quantum effects sufficiently ameliorate any chaotic classical trajectories |T7j to 
permit closed-form solutions for the wave functions, or exotic but useful forms || 
for the propagator? 

With these questions in mind, I will describe in this talkP] some new features 
of models with area-dependent potentials. As non-relativistic many-body field 
theories, I believe that area-dependent models are very interesting. I will argue 
this beginning with the classical field versions. I explicitly solve these for some 
special configurations, and then investigate general situations. The classical field 
theory problem always reduces to solving the linear Schrodinger equation in a self- 
consistent effective potential. In general, the effective potential for the V = A 
model is anisotropic and linear, while for the V = A 2 model it is anisotropic and 
quadratic. In both cases the effective potential is determined from the initial 
data by a closed set of coupled time-dependent equations. I will analyze these 
equations, especially for the case of spherically symmetric particle density, and I 
will solve them in that special situation. However, I have not completely solved 
the area-dependent models for general, anisotropic initial data. Much more work 
is needed for the V = A 2 model. Additional work is also needed to construct 
the quantum versions of these field theories. These particular models seem to 
be ripe for exploration using deformation quantization |TJ, [7|, [TI| and ideas from 
non-commutative geometry [p~5 |. 



To begin, we construct a many-body field theory on the plane with a three-body 
potential that is just the signed area of the three-body triangle (not the absolute 
value of the area). We take a collection of three types of non-identical particles, 
i.e. three particle "species", each represented by its own local field ip a , a = 1, 2, 3. 
In this case, the Hamiltonian has no lower bound for configurations with large 
but negative A, so the model is reminiscent of a linear potential single-particle 
QM. (Actually there is more than reminiscing going on here as we shall soon see.) 

the actual talk, I began with a simple model of N distinguishable point particles on 
the plane, with V = A. This is based on work in collaboration with Alexios Polychronakos 
and Cosmas Zachos We have completely understood this model's classical and quantum 

properties. Since the potential is quadratic, this is not very difficult to do, in principle, but 
nevertheless it was necessary to develop an efficient formalism to handle an arbitrary number of 
particles. We have done this using so-called cyclotomic coordinates that have previously been 
used to analyze toroidal membrane dynamics [ fl2| . I view the analysis of this V = A model as a 
precursor to point particle models with V — A 2 . These latter models are much more interesting, 
but are not so easy to solve. While I have made a modicum of progress in an unfinished attempt 
to construct a Lax pair for certain variants of V = A 2 models, in the point particle approach p3| , 
I will not discuss this here. Instead, I will shift from point particle models to field theories. 
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Explicitly, we define a multi-particle configuration by the field theory Lagrangian 



L= i f (r) (r) - # • 

a=l,2,3 ^ ( 



It should be understood that all the fields also depend on a common time, t, but 
usually we will not explicitly indicate this. We then use the obvious symmetrical 
Hamiltonian 

H = J (dr) £ ft (r) [^] ^ (r) 

a=l,2,3 

+ kfff(dr 1 )(dr 2 )(dr 3 ) A{t u t 2 ,t 3 ) |Vi (n) ^ 2 (r 2 ) ^ 3 (r 3 )| 2 
involving A, (twice) the area of the triangle formed by the three particles. 

A (ri, r 2 , r 3 ) = r x A r 2 + r 2 A r 3 + r 3 A r x . 
This may be rewritten (exactly) as 

H= J2 I (*) # W [^] ^ W + I E J(^)(5 a M+V a M-r)|^(r)| 2 

a=l,2,3 a=l,2,3 

where we have introduced the combinations 5 a [?/>] + V a [i/j] ■ r for a = 1, 2, 3 to serve 
as effective potentials for the three fields. These effective potentials are defined in 
terms of the fields as 



K M = W j £abc / / ^ ^ ^ " r *y IV* ( r 2) ^ Ml 2 

6,c=l,2,3 

SM = lke abc JJ (dr 2 )(dr 3 ) (r 2 A r 3 ) (r 2 ) ^ c (r 3 )| 2 

Note that S a and V a do not depend on ip a but do depend on the other two field 
configurations, ipb^a- Further observe that S a [tp] = if either of the two ^'s 
involved are of definite parity (either even or odd functions of r), and V* [ip] = if 
both of the two if)'s involved are of definite parity. Thus the special configurations 
where all three species have definite parity reduce to free fields! Again, this is a 
model on the plane, so each species coordinate r a is a two-vector. To construct a 
potential linear in the signed area in higher dimensions, we would need to take an 
inner product of A (ri, r 2 , r 3 ) with some constant 2-form. 

Like the point-particle theory in || , this classical field theory is always solvable, 
even in the general situation when all the fields do not have definite parity. Unlike 
|| , we work in real coordinates here. The field equations in terms of the effective 
potentials are then: 

d 1 

i— ip a (r, t) + ^^ 2 4>a (r, t) = (S a [V>] + V B [V>] • r) ip a (r, t) no sum a. 
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For each species we have just Schrodinger's equation with a linear (in r l ) poten- 
tial, albeit a non-isotropic linear potential, in general, with non-constant (in t) 
configuration- dependent coefficients. From the above definitions we obtain 

SaVl>] = \k £ e^e abc N b RlN c Ri , V* [ip] = k ^ e^e abc N b RlN c 

6,c=l,2,3 6,c=l,2,3 
M=l,2 j=l,2 



N a W a = I {dr) r j \4> a (r,t)| 2 , N a Pl = -i I {dr) ft (r,t) V'V« (r,t) 



N a = / (dr) |V a (r,t)| 



2 

j 



with no sum over a in any of these. While these are indeed configuration-dependent 
coefficients, the implicit field dependence of the coefficients is completely tractable. 
From the field equations, the coefficients in the effective potential obey simple first 
order time-derivative equations: 

^N a = 0, 4^ = —*?, 
dt dt a 2m a ' 



d Pi = -2V> = -2te« J2 z abc N b N c Rl , ^( J2 N " P < 



dt 

6,c=l,2,3 \a=l,2,3 



The last equation represents conservation of the system's total momentum, while 
the next to last equation shows that the individual particle species momenta are 
not separately conserved, in general. 

These first order equations combine to yield linear second order equations with 
constant coefficients for any given initial data (hence easily solved). 



d 



^Rl = kJ2 M MW R b > M (-)(i&) = £ij J2 z acb N c N b (no sum b) 



m dt-- 

j,b c=l,2,3 



The time- independent eigenvalues of M are and ±a/ NiN 2 N S y/ N% + N 2 + N 3 , 
with each of these three possibilities occurring twice. (Note how the three point- 
particle mechanics eigenvalues in || are obtained by setting all the iV's equal to 
one.) Therefore the i?'s, and hence all the terms in the effective potentials, can be 
solved for in terms of the initial data. The three independent classes of solutions 
for the i?'s, corresponding to the three eigenvalues of M, are functions linear in 
t, real exponentials, and oscillations. The time dependencies of the terms in the 
effective potential then follow from their expression in terms of the i?'s. Finally, 
to complete the solution, the linear Schrodinger equations for the individual fields 
must be solved using the solutions for the effective potentials. 

I believe, and assert without proof, that this last step can always be carried 
out, although there is some possibility of unusual behavior for the fields because 
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the potential terms in the Schrodinger equation are combinations of polynomials 
and complex exponentials in t. Technically, I cannot rigorously rule out some sort 
of chaotic behavior in the phases of the fields, yet, but since the behavior of all 
the field bilinears is so deterministic, including those involving spatial derivatives, 
my current opinion is that chaotic behavior in the fields, if any, must be limited to 
only time- dependent phases and is completely innocuous. (I would devote more 
time to this issue, but once again, I view this model as just a warm-up exercise 
for the A 2 model presented below.) I leave it to the reader to compare the above 
system of equations to those of the standard isotropic oscillator |^l], [2^ . 



For consistency with the treatment of the point-particle case in 0, and to warm 
the hearts of conformal field theorists, the above analysis should perhaps be re-done 
using complex coordinates. Also, a generalization to a polygon area potential, 
involving N types of particles and A (r 1; r 2 , ■ • ■ , r N ) = r x Ar 2 + r 2 Ar 3 + - • -+rjvAri, 
is straightforward. An elegant touch would be to introduce a tensor E ai "' aN on 
the particle type such that the combination A (r 1? r 2 , ■ • • , r^v) E ai "' aN is symmetric 
under any pair interchange r^, dj Tj,aj. Again the most streamlined way to 
do this on the plane is to work with complex coordinates, as was the case in the 
point-particle situation. 

We next construct an area-squared field theory, where the potential is A 2 , a 
quartic function of field coordinates. This can be incorporated into a many-body 
field theory with only a single complex field V in any number of spatial dimensions 
n. The Hamiltonian for the model is 

+ \k [ d n s x [ d n s 2 [ d n s 3 \^(s l )\ 2 \^(s 2 )\ 2 \^(s 3 )\ 2 [s 2 12 s 2 23 -(s 12 -s 23 ) 2 ] 



where s a f, = s a — Sj, for each ab pair of the identical particles' coordinate vectors. 
The equations of motion from L [if)] = j d n s i^dfif) — H[if)] are 

idfip (r) + T^V 2 ^ (r) = V [tp] (r) ip (r) 
In this case the effective potential corresponds to a non-isotropic oscillator 

V M (r)=kj d n s 2 J d n s 3 |V (s 2 )\ 2 |V (s 3 )\ 2 [s 2 l2 s 2 23 - ( Sl2 • S23 ) 2 ] | si=r 

= i^-rV + By + Z = K,ij (rV - 2R}r j ) + Z 
= K i:i (r* - B?) (r j - R j ) + Z- K^KW 

with the definitions and relations 



Kij = kj d n s 2 J d n s 3 |V (s 2 )\ 2 |V (s 3 )| (S 
= 2kN 2 ((S - R 2 ) 5ij - (S lj - RiRj)) 
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Bt = -2k J d n S 2 J d n S 3 |V (S 2 )| 2 |V (S 3 )| 2 (4 4 - S 23 S 2 ■ S23) 

= 4kN 2 (SijRj - SRi) = -ZKijRj 
Z = k J d n s 2 J d n s 3 |V (s 2 )| 2 IV (ss)| 2 (s 2 2 sl - (s 2 ■ S3) 2 ) = kN 2 (S 2 - SijSji) 

N = fd n s\^(s)\ 2 , Rl = ±fd n ss t \^(s)\ 2 , Pi [V] = ^Jd n s ($Xl>) , 
Sij = ^Id n s s iSj |V (s)| 2 , S = J2 Sjj , K i3 = 2kN 2 (SS i3 - 5 y ) . 

Note that everything in V [V] factorizes into integrals of bilinear densities. Actu- 
ally, this is not too surprising given the polynomial character and factorized form 
of the terms of the three-body potential. So the effective potential for the area 2 
model becomes 

V [V] (r) = 2kN 2 [S (r - R) 2 - (r* - IV) (r J - R j )] + 2kN 2 (S ij R i R j - SR 2 ) 
+ kN 2 (S 2 - SijSij) - 2kN 2 (R 2 r 2 - (r • R) 2 ) 

As already noted, this is a non-isotropic quadratic potential. Note the presence 
of the last term R 2 r 2 — (r • R) 2 = (r A R) 2 which is the squared area of the par- 
allelogram formed from r and R. The negative coefficient of this term and the 
potential instability it would produce for large r _L R are completely offset by the 
additional presence of 2kN 2 Sr 2 in the effective potential, along with the standard 
inequality S > R 2 . Also note that V is non-zero only through the field spreading 
out around the center of mass location. Were |V( S )| 2 a delta-function at R, all 
coefficients would vanish in V. 

Also note the absence of nontrivial "zero modes" for the quadratic part of V 
(i.e. those zero-modes that would persist even in the CM frame, and not those 
associated with the free translation of the CM). This is straightforward to see, 
even though for the field coordinates alone det (siSj) = = det (s 2 5ij — SiSj), and 
so these coordinate matrices have zero eigenvalues. However, (det S) 7^ det (S) just 
as (x 2 ) 7^ (re) 2 . In fact det (S5ij — Sy) = ^det J d n s (s 2 5ij — SiSj) |V (s)| 2 > 0. 
Moreover, in the CM det K > and all the eigenvalues of are greater than 
zero. Or equivalently, all the eigenvalues of SSij — Sij are positive. Proof: For any 
real vector v, we have Vi (SS^ — SV,) Vj = f d n s (s 2 v 2 — (s • v) 2 ) |V (s)| 2 . But 
the integrand is non-negative by the Schwarz inequality s 2 v 2 — (s ■ v) 2 > 0. Thus 
any eigenvalue A of is non- negative, A > 0. Now to rule out the case where any 
eigenvalue vanishes in the CM, we note the Schwarz inequality is only an equality 
s 2 v 2 — (s ■ v) 2 = if Si oc Vi. Therefore the only way to have a zero eigenvalue is to 
have the full integrand (s 2 t> 2 — (s ■ v ) 2 ) |V (s)\ 2 restricted to a line where Sj oc i>j. 
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Forget that! There is always some transverse spread in \ip (s)\ 2 . Any reasonable 
wave function (or classical field) in n > 2 dimensions does not have support just 
on a line. So V [ip] (r) always has a quadratic r 2 part and therefore the field is 
always localized around the CM. Or, perhaps more physically, all the particles are 
bound by the effective potential to the CM. 
For solutions the energy is 

if M = ^-NT + kN 3 (S 2 - SijSij) 
8m v ' 

using Sij = Sij — RiRj , S = Sjj , T = Tjj , with the stress-tensor defined as 

Now what is the time-dependence of the three effective potential coefficients 
Kij, Bi, and Z7 Or rather, what is the time-dependence of the quantities that 
compose them? Obviously, N, Ri, and Pi obey the equations 

d , T d „ d „ 1 „ 

—AT = , —Pi = , — Ri = — Pi 
dt dt dt 2m 

with the solutions 

N (t) = N = N , Pi (t) = P Qi = Pi, R l (t) =Roi + ^ Pit 
But what about S^l With the definition (note that Q^ ^ Qji) 

Qij [V>] = 9/ rf " SSi (^^) = ~W j (rsSi (^)-iS- 



V) 



we find 

dt Sij = 2m~ ( Qij + Qji) 

So we must include Qij into the mix of differential equations. It's time derivative 
is 

~~^Qij 2,m ^^ik-K-kj IRiBj 

where Ty [ip] is the previous stress-tensor. We are one step away from closing the 
system. We just need the final time derivative 

= ~^ (KikQkj + Kj k Q ki ) — 2PiBj — 2PjBi 

The system of equations is now closed. The problem then is to solve these equa- 
tions to obtain and hence the effective potential in terms of the initial field 
configuration . Note that only the symmetric tensor part of Qij varies with time. 
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Let's break down the tensors into somewhat more traditional combinations. 
The symmetric and antisymmetric parts of are (almost) familiar from rotations 
and (at least when traced) coordinate rescaling. 



J%j Qij Qji i -^ij Qij ~t~ Q, 



3i 

These nicely occur above with their values for (or about) the CM removed. Thus 
define 

^ij ^ ij (.RiRj ~t~ RjRi) •} 3ij (.RiRj RjRi) 

Similarly remove from the quadrupole and stress tensors the CM coordinate and 
momentum dyads. 

Sij ^ij R{Rj , Tlj ^ij ^i^j 

Incorporating momentum and angular momentum conservation, i.e. Bi = —IK^R? 
and K ik Skj — RjK ik R k = Kj k S ki — RiKj k R k (just a matrix commutator, [S, K]^ = 
0), the time evolution equations are then 

iU,. = J-iV,., 4^ = 0, 4^ = o (1) 

dt 2m dt dt 

-j-T>ij = — %j — 4:Si k K k j — ASj k K ki = — %j — 4 {S, K} ■ . 
dt m m 3 

jTij = -2K lk (V kj + J kj ) - 2K jk {V kl + J ki ) = -2 {K, V}^ - 2 [K, J\. 

where now Kij = 2kN 2 (S5ij — Sij) . Anyone skilled in the black arts of nuclear 
physics should feel comfortable with equations of this sort. 
An invariant arises from taking traces. Thus 

(tr7^ + 8mkN 2 (S 2 - trS 2 )) = . 

Now, this is essentially just the energy in the center-of-mass. 

P 2 + \x%j + 8mkN 2 (S 2 - trSfj) = trT^ + 8mkN 2 (S 2 - tr<Sj) = 8mH[iP}/N 
where 

Hty] = ^NT + kN 3 (S 2 - SijSij) . 

As a special situation, but really without loss of generality, consider co-moving with 
the system in the CM, where Ri = = Pj. Then the effective potential coefficients 
reduce to 

K i5 = 2kN 2 (S5 tJ - S^) , B t = , Z = kN 2 {S 2 - 



with the Sij, T>%^ and %j time derivatives unchanged in form. (Note in the CM 

For now, I cannot solve the non-isotropic equations ([[J). Perhaps closed-form 
solutions can be obtained in various limiting situations, but most probably the sys- 
tem is chaotic for general initial dataQ. Nevertheless, I can solve exactly a further 
specialization, which is with loss of generality. Suppose the various densities are 
actually isotropic about the CM. Then we have 

V iA = -VbiA , SiA = -S5a , Ta = -T6in , K iA = -K5a = 2kN 2 ( ] S5i 



n n n n \ n 

K = 2kN 2 {n-l)S , Z = kN 2 {^—^jS 2 

and the time derivatives are -^K^ = 2kN 2 (^-) ^P^-, etc. These may all be 
traced now w.l.o.g. to yield 

lK = -kN 2 (n-l)V, ^-Z=±kN 2 (^)sV } ^ = ^-V , 
at m at m \ n J at 2m 

and the remaining equations 

d „ 1^ , o f n — 1 \ „ d ^ , f n — 1 \ 

— V = — T - 16kN 2 )S 2 , — T = -8kN 2 SV 

at m \ n J at \ n J 

Using j- t S = the T equation becomes jfT = —8mkN 2 (^) j^S 2 which 

integrates immediately to 

T (t) = T + 8mA;iV 2 f («5 2 - S (t) 2 ) 

Inserting this into the V equation and using jrS = ^- V once again, we have a 
second-order equation for S. 

2m f-S = -% + 8kN 2 ( — ) S 2 - 24kN 2 ( — ) S 2 
at 2 m \ n J \ n J 

Thus S alone evolves as a nonlinear oscillator subject to a cubic self-interaction. 
That is, the effective potential governing the time evolution of S is V (S) = aS 3 — 
bS, where the constant, initial-data-dependent coefficients are a = 8kN 2 ( rL -^) , 
b = -^Tq + 8kN 2 Sq, and are both positive. So the potential has a local 

minimum/maximum at 3aS 2 = b or iS min / max = ±W^, hence 



1 _„ 1 In 



5 min/max - ± t / 3 5 2 + 2AmkN2 ( _ j ) % 



2 Since giving this talk, extensive numerical investigations have shown that the model is indeed 
chaotic for general initial data S . 
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Only the local minimum here is physical. Recall that we are in the CM, so that D = 
V,T = T, and most importantly the inequality S = S = J d n s s 2 \ip (s)| 2 > ,so 
we see that physically S is constrained to be positive (and S = only for \ip (s)| 2 
ultra- localized at the origin, like a delta- function ••• a collapse/condensation of 
the matter field). 

Is it possible that we encounter singularities here, in a finite time? By it's 
definition we are constrained to have S > 0, but for "special" initial data, perhaps 
the equation evolves S to zero, the boundary of the unphysical S < region, hence 
the matter field would collapse to the origin. After all, the equation for S is that 
of a nonlinear "cubic" oscillator. Do we have here a baby version of black-hole 
physics? 

The answer is no, S = is never achieved, as we now explain. The effective 
potential for S will require S > to be true for all times if and only if So > and 
£ (S) < 0, where £ is the conserved effective energy for S. 

£{S) = m(^s) +V(S)=£ 



£ = + 8kN* (*=i) 5 3 - (I T + 8kN* (*=!) S 2 ) 5 = ± (V 2 - 4 T 5 ) 

Consider the two terms on the RHS of this last expression. The standard QM 
uncertainty relation (basically the Schwarz inequality) for any two hermitean op- 
erators a and b, assuming (a) = = (b), is (a 2 ) (b 2 ) > (ab — ba)) + 

(| (ab + ba)) 2 . Each of the two terms on the RHS of the inequality are positive. 
Usually we discard the second of these RHS terms and keep the first, especially 
in the case for canonically conjugate variables with xp — px = ih. But we could 
equally well discard the first RHS term to obtain(a 2 ) (b 2 ) — {\ (ab + ba)) > 0. 
This is a strict inequality for conjugate variables. Applying this to the situation 
at hand we conclude ATS — V 2 > 0. Hence the cubic oscillator system above is 
always bound to the stable minimum, with £ < 0. This is what we wanted to 
show. 

We will discuss static solutions in more detail below. For now, we note that 
these are possible if we are at the local minimum of V (S) with T> — 0, where we 
require 

S - = S ° = ^ + ^kNl{l^r) T °- M 31 = 18™*^^)^, 

as well as T>$ = 0, of course. The fields themselves may be worked out explicitly 
in this case and are nothing but Gaussians. 
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Even in the non-static case the conserved energy for S reduces it's determination 
to quadrature and allows us to solve for the S motion in textbook fashion. 

±S = ±y/(£o -V(S))/m 

t r s ® ds i r s ® ds 



J S y/£o-V(S) V^Jso y/(S hi -S) (S-S mid ) (S-Sio) 

where for £ < 0, {Si a , S mic i, S^} are the three distinct zeroes of the cubic £ + 
bS — aS 3 with Si Q < < S mid < S (t) < S hi . The RHS here is an elliptic integral. 

r s(t) 



f yShi-Sio ds = x ( \ I Shi ~ s ™-jd \ — f ( . / s hj-s(t) I 

/o Jt{s hi -S){s-s mid ){s-s lo ) VV s hi Si„ ) \\Js hi -s mid ^y 

^ '-'mid ^ 



Shi—Slo 



K and F are the standard elliptic integrals, with F (1, z) — K (z), and F (0, z) = 0. 
Thus S (t) is an elliptic function. The solution oscillates between the distinct 
turning points S mid < S (t) < Shi with period 

u fa f Shl dS _ / 4m K ( / s hi -s mid \ 

2 t period V a V( s M -S)(S-S mid )(S-S lo ) - V IV S M Su, ) 

We may extend the analysis to higher dimensions for models with potentials 
depending on higher multi-particle coordinate forms. Assume a (d + l)-body 
potential U in n > d dimensions which is a (<i-form) 2 . That is 

F (ri, r 2 , r 3 , • • • , r d , r d+ i) = (ri - r 2 ) A (r 2 - r 3 ) A • • • A (r d - r d+ i) 



U (ri, r 2 , • • • , r d+1 ) = £;F 2 (r 1; r 2 , • • • , r d+1 ) 



L r'l«2-!d „«1 „J1 „«2 -J2 „J3 ~?3 „*d -3d 

~ ^ °hh-3d 12 1 12 1 23 1 23 1 34 1 34 ' " 1 dd+l 1 dd+1 

1 



2m 



^(r) 



if = J (dr) ij* (r) 

+ dT~[ j '" / ( rfr i) • • • ( rfr rf+i) */(ri,r 2 ,--- ,r d+1 ) |^ (r d+1 ) • • (r 2 ) ^ ( ri )| 2 

with / • • • / (dr 2 ) ■ ■ ■ (dr d+1 ) U (r, r 2 , • • • , r d+ i) |^ (r d+ i) • • -ip (r 2 )| 2 ip (r) in the 
field equation. 

If we assume spherically symmetric fields we always obtain from this isotropic 
harmonic effective potentials. 

J - J (dr 2 ) ■ ■ ■ (dr d+1 ) U (r, r 2 , • • • , r d+1 ) \ip (r d+1 ) ■ ■ -tp (r 2 )| 2 

/ M \ d ~ 1 

— k(n — i)(n — 2)---(n — d+l) (dM r 2 + M 2 ) I — - J 
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where M k = J (dr) (r 2 ) k/2 |V (r)| 2 = Q n f™ |V>| V+^dr and tt n = 2ix n > 2 /V (n/2). 
As usual the interpretation (at least in the quantum theory) is that N = M is the 
total number of particles. Thus the field equation for such spherically symmetric 
solutions becomes 

■ d , , , 1 / d 2 , , , n- 1 d , . \ Mt 1 (n- 1)! , „, , „, x , , , , 

'«* (r) + 2^ M + — W J = n I 1( „_ d) ! + *>) W 

For a Gaussian ansatz, tp(r,t) = C exp (— |cr 2 ) exp (— iEt), this field equation 
reduces to 

£ + — (c 2 r 2 - nc) = k-. — dNr 2 + 

2m v ; (n-d)! V2c/ \ 2c J 

and so we have a solution provided 

c^ 1 = A; 7 wrN d -— , E = 1 + — - 

2m (n-d)! 2*" 1 2m \ 2d 

The E above is almost, but not quite, the energy per particle. That is, the value 
of H, the total energy, is almost but not quite EN for the spherically symmetric 
solution. For the Gaussian ansatz 

ncN ( 1 
4m V d 

The above suggests what to do also in the situation where ip is not necessarily 
spherically symmetric. For any U (r, r 2 , r 3 • • • , rj, r^+i) which is a quadratic form 
in the components of r we may write quite generally 

Z + JBi + rVK tJ =/•••/ (dr 2 ) • • • (dr d+1 ) [/ (r, r 2 , • • • , r d+1 ) | V (r d+1 ) • • • V (r 2 ) | 2 

This is true for the case at hand. U (r, r 2 , r 3 • • • , r^, r^+i) is a quadratic form in the 
components of r (as well as quadratic in all the other individual r a ) even though 
it is of order 2d altogether. Hence 

K * = k ^A I- I ^ ■ ■ ■ ( rfr ^i) 44 • • • IV (rd+i) • • • V (r 2 ) | 2 

5, = ~2k5^ d J -J (dr 2 ) • • • (dr d+1 ) if r* r£ • • • r^ +1 r^ +1 | V (r d+1 ) • • • V (r 2 ) | 2 

Z = I- I ^) ■ ■ ■ ( rfr ^i) 44 44 • • ■ r^ +1 r^ +1 | V (r d+1 ) • • • V (r 2 ) | 2 

It is now straightforward to use the field equations 

(r) + ^V 2 ^ (r) = (Z + r l B t + rV'^-) V (r) 

12 



to determine a closed set of time-derivative equations obeyed by the coefficients 
in the effective potential. Once these are solved, either in special situations or 
perhaps more generally, then the field equation itself is to be solved using the time- 
dependence determined for the coefficients, in a self-consistent way. Sounds easy, 
even if it is not in practice, but perhaps the resulting non-isotropic equations can 
always be solved in closed form in some limit, such as large n. 

There was neither enough time in the talk nor enough space in this written 
version to discuss either the supersymmetric extensions of these models or their 
quantization using deformation methods |9|, ||. These subjects will be treated 
elsewhere. 

Acknowledgments 

This work was supported in part by NSF Award 0073390 and by US Department 
of Energy, Division of High Energy Physics, Contract W-31-109-ENG-38. I thank 
Cosmas Zachos and the Particle Theory Group at Argonne National Laboratory 
for their hospitality in the summer of 2001 during which a portion of this research 
was completed. 

References 



[i] 

[2] 
[3] 
[4] 

[5] 
[6] 
[7 



[9 



F Bayen, et al, Ann Phys 111 (1978) 61; ibid. 111. 

A Carnegie and I C Percival, J Phys A 17 (1984) 801-813. 

S-J Chang, Phys Rev D 29 (1984) 259-268. 

T Curtright, ICTP seminar, 1989, 

server.physics.miami.edu/~curtright/ictp89.htm" 



T Curtright, "Schrodinger's Cataplex" [|quant-ph/00 1 1 101 
T Curtright, unpublished. 



T Curtright, D Fairlie, and C Zachos, Phys Rev D58 (1998) 025002 |Eepj 
th/9711183|| . 



T Curtright and H Kocak, in preparation. 

T Curtright, A Polychronakos, and C Zachos, Phys Lett A295 (2002) 241-246 
||hep-th/01liT73 |. 



[10] T Curtright, T Uematsu, and C Zachos, J Math Phys 42 (2001) 2396-2415. 
||hep-th/001lT37| . 



[11] B de Wit, J Hoppe, and H Nicolai, Nucl Phys B305 (1988) 545-581; 
B de Wit, M Liischer, and H Nicolai, Nucl Phys B320 (1989) 135-159. 



13 



[12 

[13 

[14 
[15 

[16 



[i7; 



[19 

[20 

[21 
[22 

[23 

[24 

[25 



D Fairlie and C Zachos, Phys Lett B224 (1989) 101; 
E Floratos, Phys Lett B228 (1989) 335-340. 

S Ferrara, Christian Fronsdal, and A Zaffaroni, Nucl Phys B532 (1998) 153- 
162 [|hep-th/9802203| . 



R P Feynman, Nucl Phys B188 (1981) 479-512. 

J M Gracia-Bondfa, J C Varilly, and H Figueroa, Elements of Noncommutative 
Geometry, Birkhauser, 2001. 

G M Graf, D Hasler, and J Hoppe, "No zero energy states for the supersym- 
metric x 2 y 2 potential" [|math-ph/0109032|| ; 



J Froehlich, G M Graf, D Hasler, J Hoppe, S-T Yau, Nucl Phys B567 (2000) 
231-248 ||hep-th/9904T8H . 



M C Gutzwiller, Chaos in Classical and Quantum Mechanics, Springer- Verlag, 
1990. 



P S Howe and P C West, "Is N=4 Yang-Mills Theory Soluble?" flhep- 
th/ 96 11 074|| . 



A Jaffe and E Witten, "Quantum Yang-Mills Theory" 
| |http: / / www.claymath.org/ prizeproblems/yang_mills.pdi|| . 

S Lee, S Minwalla, M Rangamani, and N Seiberg, Adv Theor Math Phys 2 
(1998) 697-718 ||hep-th/9806074 |. 



H R Lewis, Jr., and W B Riesenfeld, J Math Phys 10 (1969) 1458-1473. 

G D Mahan, Many-Particle Physics, Kluwer Academic/Plenum Publishers 
2000. 

A M Perelomov, Integrable Systems of Classical Mechanics and Lie Algebras 
(Volume I) Birkhauser 1990. 

G K Savvidy, Nucl Phys B246 (1982) 302 and Phys Lett 130B (1983) 203; 
S G Martinyan, E B Prokhorenko, and G K Savvidy, Nucl Phys B298 (1988) 
414. 

B Simon, Ann Phys (NY) 146 (1983) 209-220. 



14 



